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We demonstrate site-resolved imaging of individual fermionic ®Li atoms in a 2D optical lattice. 
To preserve the density distribution during fluorescence imaging, we simultaneously cool the atoms 
with 3D Raman sideband cooling. This laser cooling technique, demonstrated here for the first 
time for ®Li atoms, also provides a pathway to rapid low-entropy filling of an optical lattice. We 
are able to determine the occupation of individual lattice sites with a fidelity >95%, enabling 
direct, local measurement of particle correlations in Fermi lattice systems. This ability will be 
instrumental for creating and investigating low-temperature phases of the Fermi-Hubbard model, 
including antiferromagnets and d-wave superfluidity. 


PACS numbers: 03.75.Ss, 07.60.Pb, 37.10.De 

Particle correlations reveal the underlying order of an 
interacting quantum many-body system. Strong correla¬ 
tions give rise to rich quantum many-body phenomena 
such as high-temperature superconductivity and colossal 
magneto-resistance [T]. One approach toward studying 
correlated many-body systems uses ultracold atoms to 
implement a well-understood and tunable realization of 
a particular model, and to use the behavior of the clean 
atomic system as a benchmark for theory [2]. This “syn¬ 
thetic matter” approach is especially fruitful for strongly- 
correlated fermionic systems, where, for even the simplest 
models, the sign problem of the Quantum Monte Carlo 
method precludes accurate computations of thermody¬ 
namic observables [3]. In addition to theoretical sim¬ 
plicity and tunability, ultracold atomic systems can be 
designed to have interparticle spacings of order the wave¬ 
length of visible light. By placing a quantum gas under 
an optical microscope we can therefore directly observe 
and manipulate quantum correlations at their smallest 
length scale. Such a quantum gas microscope has been 
realized for bosonic ^^Rb H ig and m atoms. In 

bosonic systems, site-resolved imaging has been used to 
study the quantum phase transition from a superfluid to 
a Mott insulator mmm and from a paramagnet to an an- 
tiferromagnet |9]. Single-site resolution also enables the 
extraction of non-local order parameters such as string 
order m and allows studies of strongly-correlated dy¬ 
namics in optical lattices [nHis]. Until very recently [H- 
IISI, however, site-resolved imaging had not been demon¬ 
strated for fermionic atoms. In Fermi-Hubbard systems, 
cold atom experiments without single-site resolution have 
observed Mott insulators [niiTH] and antiferromagnetic 
correlations da HO]. In these experiments, understand¬ 
ing of the prepared many-body state is limited by lack of 
direct access to the many-body wave function and the in¬ 
ability to locally measure correlations. The extension of 
quantum gas microscopy to fermions will provide novel 
probes for Fermi lattice systems, such as site-resolved 
spin correlation functions and local entropy measure¬ 



FIG. 1. (color online). Fluorescence image of atoms in a 
single layer of a cubic lattice obtained using Raman sideband 
cooling. The filling fraction in the center of the cloud is 40%. 
We collect approximately 750 photons per atom during a 1.9 s 
exposure. The colorbar is in arbitrary units. 


ment. 

Here, we demonstrate site-resolved imaging of 
fermionic ^Li in a 2D optical lattice with high fidelity 
[see Fig.[^. ^Li is an especially suitable species for many- 
body experiments with ultracold atoms because its light 
mass leads to fast thermalization and dynamics, and its 
broad magnetic Feshbach resonances [21] allow precise 
control of atomic interactions. The natural energy scale 
for particles of mass m, in an optical lattice with spacing 
a, is the recoil energy, Er = In? jSo ?where h is Planck’s 
constant. For many-body physics, working with a light 
atom gives an advantage because the recoil energy scales 
inversely with the mass. Experiments studying antifer- 
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FIG. 2. (color online). A schematic of the microscope. R1 and 
R2 denote our Raman beams, and OP the optical pumping 
light which co-propagates with R2. LI and L2 are addition¬ 
ally retroreflected out of the schematic to create a 3D lattice 
as described in Ref. |30| . L3 forms a lattice along z, pro¬ 
viding additional confinement during imaging. LI, L2, and 
L3 have waists of 80 /xm, 80 /xm, and 40 /xm, respectively. 
The measured point spread function, obtained by superim¬ 
posing and averaging isolated atoms, is shown in panel (b). 
The black markers are an azimuthal average of the measured 
point spread function (PSF). The red curve is the expected 
diffraction-limited Airy disk for an NA of 0.87. The inset is 
an image of the PSF. A gaussian fit to the PSF yields a full 
width at half maximum of 520 nm, compared to our lattice 
spacing of 569 nm. 


romagnetic correlations with [1^] have been limited 
by heating, owing to the intrinsic slow dynamics of cold 
atoms. The natural timescale for ^Li is 7 times faster 
than for in a system with identical lattice geometry. 
For microscopy, however, the light mass creates a chal¬ 
lenge because the recoil energy due to photon scatter¬ 
ing also scales inversely with the atomic mass, requiring 
very large trap depths for imaging. We overcome this 
challenge by implementing 3D Raman sideband cooling 
[ 2 M 29 ] for ^Li atoms in a 2.4 mK deep optical lattice. 

Atoms are trapped in a vacuum glass cell, 9.9 /xm be¬ 
neath the surface of a superpolished substrate, in the 
object plane of a 0.87 numerical aperture (NA) imag¬ 
ing system. Our imaging system combines a long work¬ 
ing distance microscope objective (Optem 20X, NA=0.6) 
with a hemispherical lens to enhance the NA. We com¬ 
pensate spherical aberration with a phase plate in the 
imaging system. We image the atomic fluorescence onto 
the photocathode of a gateable intensified CCD camera 


(Andor iStar 334T) with a magnification of 170. We 
achieve diffraction-limited resolution, shown in Fig.|^b). 
The full width at half maximum from a Gaussian fit to 
the measured point spread function is 520 nm compared 
to a lattice spacing of 569 nm. 

Atoms in an equal mixture of \F = 1 / 2 , mp = iV^) 
the electronic ground state are loaded from a single layer 
of a ID “accordion lattice” with tunable spacing into a 2D 
optical lattice (see Ref. [30]). Lattice beams LI and L2 
[see Fig. [^a)] form radial lattices along x and y, respec¬ 
tively, with 569 nm spacing. LI and L2 also each form 
an axial lattice along i with 1.48 /xm spacing. During 
the initial lattice loading, LI and L2 are each ramped 
up in 100 ms to give radial lattice depths of 30 Er^rad, 
where tunneling is suppressed. For imaging, we intro¬ 
duce an additional lattice along z, with 534 nm spacing, 
formed by L3. All lattices are derived from 1064 nm 
light. Just before imaging, LI, L2, and L3 are ramped in 
100 ms to give nearly-degenerate on-site trap frequencies 
of = 27 r X (1.25 MHz, 1.25 MHz, 1.47 MHz), 

calibrated using lattice modulation spectroscopy. 

To keep the atoms pinned to their lattice sites during 
fluorescence imaging we must simultaneously cool them. 
Previous quantum gas microscopes have used a polar- 


(®) Raman Pump (b) 



FIG. 3. (color online). Pulsed Raman sideband imaging. A 
Raman transition drives atoms into the |2^Si/2(E = 2 / 2 )) hy- 
perfine manifold, removing one vibrational excitation. Atoms 
are then optically pumped back into the | 2 ^Si/ 2 (E = V^)) 
manifold while simultaneously switching on the intensifier of 
an intensified GGD camera to collect the photons scattered 
during pumping. A spectrum, taken by driving a Raman 
transition with a 200 /xs long pulse and then imaging the 
| 2 ^Si/ 2 (E = Y2,?7Xf = “V^)) state is shown in panel (b), with 
the red line denoting the two-photon detuning during imag¬ 
ing. The timing of two imaging pulses is shown in panel (c). 
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ization gradient cooling (PGC) scheme for imaging 
nig. PGC is not suitable for sub-Doppler cooling of ^Li 
due to the unresolved hyperfine splitting in the excited 
state m- Sisyphus cooling has been demonstrated for 
^Li in free space [32] and gray-molasses cooling has been 
demonstrated for ^Li both in free space and in an optical 
dipole trap [33|. These cooling techniques, however, have 
not yet been extended to the tightly-confined regime of 
optical lattices with ^Li. We use Raman sideband cooling 
because it does not rely on resolved hyperfine structure 
and has been demonstrated to cool a variety of atomic 
species to the motional ground state in optical lattices 
[231 El, optical tweezers nanzi, and ion traps [22], as 
well as to image ^^Rb atoms in optical tweezers [28] and 
optical lattices [29]. 

To image the atoms we collect the photons scattered 
during optical pumping in the pulsed Raman sideband 
cooling scheme shown in Fig. The imaging is per¬ 
formed at a magnetic field of < 20 mG. First, a Raman 
transition drives the atoms into | 2 ^Si/ 2 (T = 3 / 2 )), remov¬ 
ing one vibrational excitation. The Rabi frequency for a 
Raman cooling transition on the lowest motional side¬ 
band for lattice axis u is given by 77zyOc, where 77 ^^ = 
6kjyXjy = (0.47,0.47,0.15). Here, is the Lamb-Dicke 
parameter for the Raman transition, 6ky is the projection 
of the difference in the Raman beam wave vectors along 
the lattice axes, is the harmonic oscillator length, and 
lie = 27r X 160 kHz is the two-photon Rabi frequency 
on the carrier. The Raman beams have linear polariza¬ 
tion to avoid effective magnetic fields. During the Ra¬ 
man pulse, the camera intensifier is gated off to suppress 
background from the Raman light. After a 5 /is Raman 
pulse, the atoms are pumped with resonant light through 
| 22 Pi/ 2 (F = 3/2)) back into the \ 2 ^Si/ 2 {F = y 2 )) dark 
state at a rate of ^1.5 x 10^ s“^ for 20 /is, complet¬ 
ing one imaging pulse. The camera intensifier is gated 
on during the optical pumping step to collect the scat¬ 
tered photons and form an image. To obtain one image 
with ^750 photons collected per atom, we apply 6.4 x 10^ 
imaging pulses over 1.9 s. 

For efficient cooling the system must be in the Lamb- 
Dicke regime, 770 P = kopXu 1, where the optical 
pumping process preserves the vibrational state with high 
probability. Here, ?7 op ~ 0-31 is the Lamb-Dicke parame¬ 
ter for the pumping process, and kop is the magnitude of 
the wavevector for the pump light. Achieving the Lamb- 
Dicke regime for ^Li requires MHz-level trap frequencies, 
which are atypically large for neutral atom experiments 
m- The small lattice beam waists in the experiment 
cause inhomogeneity of the trap frequency over the sam¬ 
ple size. The lattice along z has the largest inhomogene¬ 
ity, with the trap frequency varying by 120 kHz over a 
radius of 30 lattice sites. We have found that the imag¬ 
ing works optimally for Raman pulse durations of 5 /is, 
where Fourier-broadening exceeds the inhomogeneity in 
trap frequency. Additionally, we find a strong depen¬ 


dence of the imaging fidelity on the detuning of the op¬ 
tical pumping light [see Fig. |^c)]. The optimal pump 
detuning is in agreement with the expected shift of the 
pump resonance due to the AC stark shift in the lattice, 
based on the polarizabilities calculated in |34| . 

We reconstruct the atom distribution in the lattice by 
fitting images to a lattice of the measured point spread 
function (PSF) [see Fig. Qa)]. The PSF of the imag¬ 
ing system and the lattice geometry are determined once 
from images of a sparsely-filled lattice and used for fitting 
subsequent images. Each image is divided into 10 x 10- 
site subregions, and each subregion is fitted with PSF am¬ 
plitudes for each site, a uniform background, and a global 
2D coordinate offset as fit parameters. A threshold is 
then applied to the fitted amplitudes to determine which 
sites were occupied. A histogram of the fitted amplitudes 
[see Fig. ib)] shows a bimodal distribution with the 
peaks corresponding to unoccupied and occupied sites. 


(a) 





Raman detuning (MHz) 


Fit amplitude (a.u.) 



Number of pulses (10"^) 


FIG. 4. (color online). Site-resolved imaging with high fi¬ 
delity. Images obtained by applying 6.4 x 10^ Raman imag¬ 
ing pulses are fit to a lattice of the measured point spread 
function, with the amplitudes for each lattice site as fit pa¬ 
rameters (a). Whether a site is occupied is determined by 
applying a threshold to the fitted amplitudes, with occupied 
sites denoted by white dots. An averaged histogram of fitted 
amplitudes is shown in panel (b). We optimize the imaging by 
taking two subsequent frames, with additional imaging pulses 
in between, and looking at the fraction of atoms that remain 
pinned to their sites between the two images (c, d). By vary¬ 
ing the number of pulses between frames and applying a linear 
fit to the pinned, hopping and lost fractions, we extract loss 
and hopping rates (e). As expected, within the statistical un¬ 
certainty of the fit, the x-intercepts of the hopping and lost 
fractions occur at 6.4 x 10^ pulses (equivalent to one frame). 
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We do not observe peaks corresponding to more than one 
atom per site because pairs of atoms are ejected during 
imaging due to light-assisted collisions [4]. Both the mo¬ 
tion of atoms between lattice sites during imaging and 
the quality of the image fit contribute to the imaging 
fidelity. By simulating images—taking into account pho¬ 
ton shot noise, camera noise, image background, and the 
measured variance in atom fluorescence—we evaluate the 
accuracy of the density reconstrunction algorithm alone, 
isolated from the effects of atomic motion. The accuracy 
is determined by comparing the known density distribu¬ 
tion in simulated images with the results from applying 
our fitting algorithm to the same images. For a lattice 
with 20% of the sites occupied, we find that the algorithm 
correctly identifies occupied sites (98.7±0.5)% of the time 
and correctly identifies unoccupied sites (99.7 ±0.2)% of 
the time. 

To study atom hopping and loss due to the imag¬ 
ing, we take two images with 6.4 x 10^ Raman imag¬ 
ing pulses each and apply a varying number of Raman 
imaging pulses applied in between them. By comparing 
the reconstructed atom distribution of the two frames, 
we determine the fraction of atoms that stay pinned to 
their sites, hop between sites, and are lost from the im¬ 
age (/p, /h, and /i). Loss can be caused by atoms leav¬ 
ing the region of analysis, hopping along i, or leaving 
the trap. The Raman imaging parameters are optimized 
on the pinned fraction [see Fig. [^c, d)]. Fig. Qe) shows 
/p, /h, and f\ versus the total number of imaging pulses 
for optimized imaging parameters. By applying a lin¬ 
ear fit to these data, we determine rates that we use to 
get the expected pinned, hopping, and lost fractions for 
a single image with 6.4 x 10^ pulses. For a single im¬ 
age we have /p = (95.1 ± 1.2)%, /h = (2.3 ± 1.3)%, and 
/i = (2.6 ± 1.7)% in a 20x20-site analysis region. A 
negative lost fraction corresponds to atoms entering the 
region of analysis. In a lattice with unity filling, each 
hopping event will cause the loss of two atoms due to 
light-assisted collisions on doubly occupied sites. Atoms 
have uniform probability of hopping at any time during 
the imaging process. From the histogram, we see that 
an atom which hops in the last half of the imaging se¬ 
quence will still be counted by the density reconstruction 
algorithm. We estimate the probability of accurately de¬ 
termining the occupation of a lattice site to be >95%. 

In conclusion, we have demonstrated site-resolved de¬ 
tection of fermionic ^Li in a 2D optical lattice with high 
fidelity using 3D Raman sideband cooling. The micro¬ 
scope will provide exquisite control of optical potentials, 
enabling single-atom addressability and perhaps a route 
to lower entropy samples [35l [36] . The extension of quan¬ 
tum gas microscopy to fermionic systems enables local 
measurement of particle correlations and will allow new 
experimental comparisons to the predictions of interact¬ 
ing quantum many-body models. 
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